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ABSTRACT 

Linearized  cavity  flow  theory,.  In  Gonjune— 
tlon  with  a  mapping  technique.  Is  used  to  develop 
general  expressions  for  the  characteristics  of  super- 
cavltating  or  fully  ventilated,  two  dimensional  hydro¬ 
foils  with  prescribed  pressure  distributions,  designed 
for  operation  near  a  free  surface „  It  Is  assumed 
that  the  hydrofoils  are  operating  at  zero  cavitation 
niunber  and  that  the  Proude  number  -  based  on  depth  - 
Is  very  large.  The  general  expressions  are  used  to 
derive  the  lift,  cavity  drag  and  shape  of  hydrofoils 
composed  of  2-,  3-j  and  5-term  and  constant  pressure 
camber  configurations,  combined  with  angle  of  attack 
and  qiiasl-parabollc  thickness  and  designed  for  opera¬ 
tion  at  specific  depths.  Examples  of  numerical  re¬ 
sults  are  given  and  the  effect  of  foil  strength  Is 


discussed  and  evaluated. 
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NOTATION 


k,  a 


Constants  In  the  conformal  transformation,  depend¬ 
ing  on  depth 


Drag  Coefficient  = 


W  0 


Lift  Coefficient  = 


Moment  Coefficient  = 


M 


X  tt2  2 

tP  U  c 


c 

Cj5(h,k^,6,T) 


C^  (h,  j  6 ) 


Cj^(h,kj^,6) 


Chord  length 

Drag  Coefficient  of  a  foil  designed  for  depth- 
chord  ratio  h,  with  camber  distribution  type  1, 
camber  index  k^^,  design  angle  of  attack  6  and 

qua si -parabolic  thickness  coefficient  t 

Lift  Coefficient  of  a  foil  designed  for  depth- 
chord  ratio  h,  with  camber  distribution  type  1, 
camber  index  k^^,  and  design  angle  of  attack  5 

Moment  Coefficient  of  a  foil  designed  for  depth- 
chord  ratio  h,  with  camber  distribution  type  1, 
camber  index  k^,  and  design  angle  of  attack  6 


D 

h 


L 


Two  dimensional  cavity  drag 
Depth  of  submersion  /  chord 

Camber  index  of  camber  distribution  type  1.  1 

takes  on  the  values  1,  2,  3  or  5  to  denote  that 
the  camber  distribution  Is  of  the  constant  pres¬ 
sure,  two-term,  three-term  or  five-term  type,  re¬ 
spectively 

Two  dimensional  lift  on  the  hydrofoil 


M  Two  dimensional  moment  about  the  leading  edge  of 

the  hydrofoil,  positive  In  the  direction  tending  to 
pitch  the  nose  up 
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p  Local  pressure  on  the  bottom  surface  of  the  hydrofoil 

p  Cavity  pressure 

c 

static  pressure  In  the  free  stream 

Ap  Pressure  difference  between  the  upper  and  lower  surfaces 

of  the  hydrofoil 

5  Leading  edge  suction  force  on  a  fully  wetted  airfoil, 
positive  In  the  direction  of  U 

OO 

u  The  X  component  of  the  perturbation  velocity  In  the 

hydrofoil  plane,  expressed  as  a  fraction  of 

u  The  X  component  of  the  perttirbation  velocity  in  the  air¬ 

foil  plane,  expressed  as  a  fraction  of 

V  The  y  component  of  the  perturbation  velocity  In  the 
hydrofoil  plane,  expressed  as  a  fraction  of 

V  The  y  component  of  the  perturbation  velocity  in  the 
airfoil  plane,  expressed  as  a  fraction  of 

X  A  space  coordinate  In  the  hydrofoil  plane,  parallel  to 

Its  origin  coincides  with  the  leading  edge  of  the 

foil.  Distance  along  the  chord  of  the  hydrofoil  ex¬ 
pressed  as  a  fraction  of  chord  length 

y  Vertical  space  coordinate  In  the  hydrofoil  plane, 

expressed  as  a  fraction  of  the  chord  length 

y  Ordinate  of  the  bottom  surface  of  the  hydrofoil  expressed 

as  a  fraction  of  the  chord  length 

y^  Ordinate  of  the  top  cavity  boundary,  expressed  as  a  frac¬ 

tion  of  the  chord  length 

Z  Hydrofoil  section  modulus 

—  -5 

Z  Non-dimensional  section  modulus  =  Z/c 

z  A  complex  variable  =  x  +  1  y 

y  Vortex  strength 

6  Design  angle  of  attack  of  the  hydrofoil 

^  A  complex  variable  =  ^  +  1 
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Tj  Imaginary  coordinate  In  the  airfoil  plane 

6  A  dummy  variable  in  the  airfoil  plane,  defined  by 

4  =  i(l-coa0) 

4  Real  coordinate  in  the  airfoil  plane 

i'  A  dummy  variable  along  the  ^  axis 

p  Mass  density  of  incompressible  fluid 

a  Cavitation  number  =  (p  -  P  )  /  ipU 

00  ^ 

T  Quasi-parabolic  thickness  coefficient 

Tlie  quantities  A,  B,  C,  D,  E,  P,  G,  H,  J,  K,  L,  M, 
g',  h',  J',  K',  l',  m',  R,  Q  and  b  are  defined  on  pages  l8  -19 
Lg  Dllogarlthm,  defined  on  page  23 

The  quantities  P  and  Q  are  defined  on  page  22 
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1. 

2. 
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4. 

5. 

6. 


7. 


8. 

9. 


(a) 

(b) 


Supercavi bating  or  ventilated  flow  near  a  free 
aurface  at  a  =  0. 

The  linearized  version  of  a  supercavitating  hydro¬ 
foil  near  a  free  surface  transformed  into  its 
equivalent  airfoil;  a  =  0. 

Transformation  parameters. 

The  X  -  4  transformation. 

Variation  of  the  optimian  lift-drag  ratio  with  depth. 

Pressure  distribution  on  the  bottom  of  the  equiva¬ 
lent  airfoil  for  fovir  types  of  camber  distribution. 

The  Influence  of  depth  on  the  lift  coefficient, 
drag  coefficient  and  the  location  of  the  center  of 

pressure  i’pr~)  due  to  camber,  angle  of  attack  and 

quasi -parabolic  thickness  (x). 

Lift-drag  ratio  of  the  two-term  camber  without 
angle  of  attack  and  quasi-parabolic  thickness. 

Bottom  surface,  y^,  and  top  cavity  boundary,  y^, 
of  four  camber  configurations  designed  for  specific 
depths . 

One  chord  depth. 

Infinite  depth. 

Bottom  sxu’face  y  ,  and  top  cavity  boundary  y  , 

O  0 

of  quasi-parabolic  thickness,  at  various  depths 
of  submersion. 
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Plgure  Number 

11.  Top  cavity  boundary  due  to  an  inclined  flat  plate 

at  various  depths. 

12»  Llft-^Drag  ratio  versua  lifi7  cnoefflclent  or  a  twO- 

term  foil  with  angle  of  attack  and  quasi-parabolic 
thickness,  t  =  0.004;  h  =  1. 

13"  Typical  two-dimensional,  two-term  hydrofoil  de¬ 

signed  for  operation  at  a  depth  of  one  chord. 
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THE  LINEARIZED  THEORY  FOR  SUPERCAVITATINQ  HYDROFOILS 
OPERATING  AT  HIGH  SPKBPS  NEAR  A  FREE  SURFACE 

INTRODUCTION 

Increasing  Interest  in  the  design  of  high  speed  hydrofoil 
craft  has  made  it  necessary  to  review  some  of  the  basic  contribu¬ 
tions  in  the  field  of  supercavltating  or  ventilated  lifting  sur- _ 

faces  and  to  investigate  the  effect  of  the  proximity  of  a  free 
surface  on  their  design.  The  exact  solution  of  the  problem  of  a 
supercavltating  inclined  flat  plate  near  a  free  surface  was  ob- 
talned  by  Green  (l)  and  was  necessarily  given  in  terms  of  the  spray 
thickness  rather  than  the  depth  of  submersion.  This  solution  was 
discussed  by  Johnson  (2)  who  presented  it  in  practicable  form  to¬ 
gether  with  experimental  relationships  between  the  spray  thickness 
and  the  depth  of  submersion.  It  is  evident,  hovjever,  that  in  view 
of  the  large  cavitation  drag  on  a  flat  plate  at  both  finite  and 
infinite  depths,  it  is  of  great  importance  to  specify  low  drag 
(i.e.,  cambered)  hydrofoil  sections. 

1‘ulln  and  Burkart  (3)  first  showed  that  the  problem  of  the 
supercavltating  hydrofoil  could  be  converted  into  the  problem  of  a 
thin  airfoil  whose  lower  surface  pressures  are  positive  or  zero  and 
that,  for  reasons  of  hydrodynamic  efficiency,  the  center  of  pres¬ 
sure  of  the  foil  should  be  as  close  to  the  trailing  edge  as  possible. 
They  Introduced  a  low  drag  supercavltating  camber  configuration  by 
specifying  the  pressure  distribution  on  the  equivalent  airfoil  as 
the  first  two  terms  of  a  trigonometric  series  and  optimizing  its  co¬ 
efficients,  and  proceeded  to  calculate  the  characteristics  of  specific 


* 

Numbers  in  parenthesis  refer  to  the  list  of  References. 
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hydrofolls,  designed  for  operation  at  infinite  depth,  by  combining 
that  camber  configuration  with  angle  of  attack.  Anticipating  a 
possible  Increase  in  hydrodynainlc  efficiency,  Johnson  (4)  intro¬ 
duced  camber  configurations  derived  from  pressure  distributions  de¬ 
fined  by  the  first  three  and  the  first  five  terms  of  a  trigonometric 
series  whose  coefficients  are  also  optimized  so  as  to  give  zero  or 
positive  pressures  on  the  lower  surface  and  a  center  of  pressure 
as  far  aft  as  possible.  Ifydrofolls  incorporating  these  types  of 
camber  configuration  have  since  come  to  be  known  as  Two-,  Three- 
and  Five-Term  Foils  and  are  very  useful  in  infinite  depth  applica¬ 
tions,.  Specifically,  two-term  foils  are  being  used  in  the  design 
of  superca vita ting  propellers  (5). 

If  a  hydrofoil  designed  for  infinite  depth  is  to  be  operated 
at  a  finite  depth  the  necessity  for  maintaining  the  condition  of 
shock  free  entry  (l.e.  the  dividing  streamline  coincides  with  the 
leading  edge)  requires  that  the  foil  be  operated  at  an  angle  of 
attack  which  is  smaller  than  the  original  (infinite  depth)  design 
angle.  The  lift  and  the  chordwlse  pressure  distribution  due  to 
camber  also  change  as  the  depth  of  submersion  is  varied.  An  approxi¬ 
mate  method  for  calculating  the  effect  of  depth  on  angle  of  attack 
and  on  the  lift  due  to  camber  has  been  given  by  Johnson  (6)  who 
also  performed  experiments  which  corroborated  the  theoretical  pre¬ 
dictions  almost  completely. 

Johnson  resorted  to  an  approxlmfitlon  because  of  the  very 
complicated  Integrals  which  arise  in  applying  even  the  linear  theory 
to  the  problem  of  determining  the  finite  depth  characteristics  of 
a  superca vl bating  hydrofoil  of  specified  shape.  The  results  he 
obtained  are  adequate  and  very  useful  in  ascertaining  the  effect 
of  depth  variation  on  the  performance  of  a  given  foil.  It  is  pos¬ 
sible,  however,  to  design  hydrofoils  for  specified,  finite  depths 
without  having  to  resort  to  approximations  to  the  linear  theory. 
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Thls  Is  done  by  specifying  the  pressure  distribution  that  la  to  re¬ 
sult  at  the  eslred  operating  depth,  rather  than  specifying  the 
shape  of  the  hydrofoil.  This  approach  results  In  hydrofoils  which 
are  designed  for  specific  depths  and  deliver  hydrodynamic  perform¬ 
ance  that  Is  somewhat  superior  to  that  of  foils  designed  for  In¬ 
finite  depth  but  operating  at  finite  depth.  In  the  present  work 
the  linearized  theory  (3)  Is  used  to  derive  general  expressions.  In 
terms  of  an  arbitrary  pressure  distribution,  for  the  shape  and  per¬ 
formance  characteristics  of  any  supercavltatlng  (0=0)  hydrofoil 
designed  for  any  specified  depth.  In  addition  to  the  limitations 
of  linear  theory  the  only  conditions  atta-ched  to  the  validity  o-f 
these  expressions  Is  that  the  press\ires  be  positive  over  the  entire 
bottom  surface  of  the  foil,  that  the  cavity  begin  at  the  leading 
edge  and  that  the  cavitation  number  based  on  cavity  pressure  be 
zero. 

The  general  expressions  are  used  to  calculate  the  lower  sur¬ 
face  shape,  cavity  shape,  lift  and  drag  of  two-,  three-  and  five- 
term  hydrofoils  and  of  a  hydrofoil  on  which  the  pressure  distribu¬ 
tion  is  constant  over  the  chord.  Prom  the  results  of  these  calcu- 
latlons  and  from  a  consideration  of  a  hypothetical  foil  on  which 
the  pressure  distribution  is  that  due  to  a  single  vortex  It  Is 
evident  that  the  hydrodynamic  efficiency  of  a  foil  depends  to  a 
very  great  extent  on  the  location  of  Its  center  of  pressure.  Prom 
the  point  of  view  of  Increasing  the  lift-drag  ratio  It  would  be  de¬ 
sirable  to  have  the  center  of  pressure  as  close  as  possible  to  the 
trailing  edge.  Strength  considerations,  however.  Indicate  that 
the  center  of  pressure  should  be  located  near  the  foil  leading  edge 
and  dictate  that  the  loading  at  the  nose  not  vanish.  The  two-, 
three-  and  five-term  camber  configurations.  If  used  by  themselves 
without  additional  angle  of  attack,  would  deliver  very  high  lift- 
drag  ratios  but  would  leave  much  to  be  desired  from  the  point  of 
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vlew  of  strength.  Specifically,  the  value  of  L/D  for  a  given  lift 
Increases  as  the  center  of  pressure  Is  moved  to  the  rear,  but  In 
all  three  camber  shapes,  when  the  lift  Is  produced  by  camber  alone, 
the  location  of  the  upper  cavity  boundary  near  the  leading  edge  Is, 
according  to  the  linear  theory,  below  the  theoretical  bottom  surface 
of  the  foil.  Tlae  constant  pressure  distribution  results  In  camber 
configurations  whose  strength  Is  somewhat  Improved  but  whose  lift- 
drag  ratios  are  relatively  poor. 

The  first  remedy  that  suggests  Itself  Is  the  standard  one  of 
combining  the  various  camber  distributions  with  angle  of  attack. 

The  high  leading  edge  pressures  associated  with  angle  of  attack 
have  the  effect  —  at  the  cost  of  reduced  lift-drag  ratios  —  of 
greatly  enlarging  the  cavity  thickness  and  hence  (since  the  foil  can 
be  thickened  to  fill  the  entire  cavity  without  affecting  Its  hydro- 
dynamic  characteristics)  result  In  increased  strength.  The  lin¬ 
earized  characteristics  of  a  flat  plate  operating  near  a  free  sur¬ 
face  are  calculated  by  using  Munk's  Inversion  formula  (7)  to  derive 
the  pressure  distribution  and  Incorporating  it  into  the  above  men¬ 
tioned  general  expressions  for  lift,  drag  and  upper  cavity  shape. 

In  addition,  the  strength  characteristics  can  be  Improved  -- 
at  moderate  cost  in  drag  —  by  superimposing  upon  the  foil  the  pres¬ 
sure  distribution  generated  by  a  single  vortex  located  at  the  lead¬ 
ing  edge  of  the  equivalent  airfoil.  The  bottom  shape  associated 
with  this  leading-edge-concentrated  pressure  Is  a  semi -parabola  in 
foils  designed  for  Infinite  depth  and  a  similar,  quasl-parabollc 
configuration  In  foils  designed  for  finite  depths.  Thus  the  pres- 
s\ire  concentration  at  the  nose  has  the  effect  of  rounding  the 
otherwise  very  sharp  leading  edge  In  addition  to  Increasing  the 
overall  area  of  the  foil  section. 

It  will  be  shown  how  the  performance  characteristics  and 
shapes  due  to  camber  (e.g.,  two-term,  three- term,  five-term  or 
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constant  pressure)  are  to  be  combined  with  those  due  to  angle  of  at¬ 
tack  and  quasl-parabollc  thickness  to  obtain  the  shape  and  charac¬ 
teristics  of  a  foil  composed  of  a  set  of  these  three  parameters. 

As  In  the  work  of  Tulln-Burkart  (3)  and  Johnson  the 

method  selected  for  the  treatment  of  the  hydrofoil  problem  Is  the  use 
of  the  concept  of  an  equivalent  airfoil  derived  from  a  suitable  con¬ 
formal  transformation. 


THE  HYDROFOIL-AIRFOIL  EQUIVALENCE 


In  considering  the  problem  of  a  supercavltatlng  or  fully  ven¬ 
tilated  hydrofoil  at  a  finite  deptlx-helavj  a  free  surface  it  Is  assiimed 
that  (a)  the  limitations  of  linear  theory  apply  (i.e.,  the  velocity 
perturbations  are  small  in  comparison  to  the  forward  speed  of  the 
foil),  and  (b)  the  cavitation  number  a  is  zero.  Asumptions  (a)  and 
(b)  are  the  same  as  those  used  In  the  design  of  Infinite  depth  hydro¬ 
foils.  In  addition  it  is  postulated  that  (c)  the  Proude  number  based 
on  depth  Is  very  large  so  that  on  the  free  sxirface  the  velocity  per¬ 
turbation  In  the  X  direction  (the  direction  of  the  stream)  is  zero. 

With  these  assumptions,  the  problem  of  the  thin  hydrofoil  near 
a  free  surface,  schematically  illustrated  in  Figure  1,  can  be  repre¬ 
sented  by  a  seml-lnflnlte  slit  on  the  x  axis  in  the  complex  z  =  x  +  iy 
plane,  the  free  surface  being  represented  by  the  line  y  =  h.  This 
configuration,  shown  in  Figure  2  together  with  the  relevant  boundary 

conditions,  is  a  simply  connected  polygon  B  ,  C  ,  C,  D,  E,  B  which 

00  00  00  00 


lends  Itself  to  a  Schwarz  Chrlstoffel  transformation  of  the  type 
dC  C+a 

where  z  =  x  +  ly  Is  the  complex  space  coordinate  in  the  hydrofoil 

plane, 

=  I  +  iq  Is  the  complex  space  coordinate  In  the  airfoil 
plane,  and  . 


[1] 
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A  and  a  are  constants  depending  on  the  depth/chord  ratio  h. 
This  transformation  maps  the  entire  z  plane  into  the  lower  half  of 
the  ^  plane. 

With  the  above  mentioned  boundary  conditions  and  the  con¬ 
straints  that 

when  X  =  0  ^  =  0 

and  when  x  =  1  ^  =  1  , 

the  Integration  of  the  transformation  [l]  and  Its  separation  Into 
the  real  and  Imaginary  parts  yield  the  result 

X  +  10  =  aA[  ^  -  ln(l  +  ]  where  ^  <  0 

^  [2] 

X  -  10  =  aA[  -  ln(l  +  ■^)  ]  where  f  >  0 
a  a 


v;here 


and 


aA-r  =  h 


In 


H-a 

a 


[3a] 

[3b] 


Plots  of  a  and  A  vs.  depth-chord  ratio  and  of  —  vs.  ^  are  shown 

aA  a 

in  Figures  3  and  4, 

As  shown  in  Figure  2  and  indicated  in  equation  [2],  the  upper 
half  of  the  z  plane  maps  into  the  negative  ^  axis  and  the  lower 
half  of  the  z  plane  maps  into  the  positive  ^  axis.  Thus  the  abclssa 
of  a  point  (x,  y^)  on  the  upper  cavity  streamline  in  the  z  plane  Is 
represented  by  a  negative  4  in  the  C  plane  and  the  abclssa  of  a  point 
(x,  y^)  on  the  lower  cavity  streamline  (l.e.,  the  lower  surface  of 
the  foil  and  the  lower  boundary  of  the  trailing  cavity)  Is  represented 
by  a  positive  ^  in  the  ^  plane.  Equations  [3al  and  [3b]  together 


» 

All  distances  are  normalized  with  respect  to  the  chord  length  of  the 
hydrofoil  and  all  perturbation  velocities  are  normalized  with  respect 
to  the  speed  of  the  stream  at  Infinity,  U  . 
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constitute  an  Implicit  solution  for  a  and  A  in  terms  of  the  depth- 
chord  ratio  h. 

At  the  limits  of  zero  and  infinite  depth  the  asymptotic 
values  of  the  constants  and  the  transformation  are: 


"wheTi  h  -+  0 


a  -*>  0 
A  1 

x-iO  -*•  ^  where  ^  >  0 


Note  that  on  the  surface  (h  =  O)  there  is  no  distortion  between  the 
airfoil  and  the  hydrofoil. 


when  h  -►  00  ^  x+10  -*■  where  ^  <  0 

^  x-iO  -*>  4^  where  ^  >  0 

Thus,  at  infinite  depth  the  transformation  [2]  reduces  to 
the  one  used  by  Tulin  and  Burkart  (3). 

LIFT,  DRAG,  MOMENT  AND  CAVITY  SHAPE 

The  lift  on  the  hydrofoil  and  the  moment  about  lus  leading 
edge  may  be  written  as 


Apdx 


Apxdx 


(positive  in  the  direction  tending  to  pitch  the  nose  up) 
where  Ap  =  the  pressure  difference  between  the  face  and  back  of  the 
Hydrofoil. 
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Using  the  linearized  Bernoulli  equation 
P  +  =  Pco 

where 

.  p  =  the  mass  density  of  the  fluid 

=  the  velocity  very  far  upstream  and  downstream  of 
the  foil 

u  =  u(x)  =  the  ratio  of  the  perturbation  velocity  on 
the  lower  (pressiire)  side  of  the  hydrofoil  to 
the  free  stream  velocity 

and  p  =  the  pressure  on  the  lower  (pressure)  side  of  the 

foil, 

and  noting  that  since  the  cavitation  number  is  zero  the  pressure 
on  the  upper  siu?face  of  a  superca  vita  ting  foil  is  p^,  we  get 

Ap  =  P-P„  =  -  p  uf u  . 

In  the  linearized  cavity  flow  theory  lines  of  constant  per¬ 
turbation  velocities  remain  unchanged  when  transformed  from  the  z 
to  the  ^  plane  (and  vice  versa).  Consequently  the  conformal  trans 
formation  has  no  effect  on  the  perturbation  velocities  u  and  v, 
other  than  to  distort  the  x  location  at  which  a  given  velocity 
exists  (l.e.j  the  pressure — proportional  to  u--at  a  point  x  on  the 
hydrofoil  is  exactly  the  same  as  the  pressure  at  the  corresponding 
point  ^  on  the  airfoil)^  so  that 
u(x)  =  u  (^) 

(unbarred  quantities  refer  to  the  hydrofoil;  barred  quantities,  to 
its  equivalent  airfoil). 

Combining  the  real  part  of  equation  [l]  and  equations  [2], 

1  2 

[9]^  and  [10],  and  normalizing  with  respect  to  we  get,  for 

00 
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the  lift  and  moment  coefficients. 

Cl  = -2  u(Od|  [111 

and 

^  ^ 


In  the  limits  of  h  =  0  and  h  -  <»  these  expressions  reduce  to 


-2 

0 

u(?)dC 

(h=0) 

[11a] 

=L  = 

-2 

0 

2iuii)di 

(h=oo) 

[lib] 

and 

2 

0 

^u(4)d? 

(h=0) 

-*• 

[12a] 

2 

2ih[(i)d^ 

(h^^ ) 

[12b] 

0 

In  linear  theory,  the  drag 

on  an  element  of  chord  of  a 

hydrofoil  may 

be  expressed  as  the 

product  of  the  lift  on  and  the 

slope  of  that  element.  Using  the  linearized  boundary  condi¬ 

tions  (shown  in  Figure  2)  and  the  fact  that  angles  (e.g.,  are 
not  affected  by  the  trans forma tion,we  may  write 

g  =  v(x)  =  v(C)  [13] 

where  v(x)U^  =  the  perturbation  velocity  in  the  y  direction  on  the 
bottom  of  the  hydrofoil, 

and  v(4)U^  =  the  perturbation  velocity  in  the  tj  direction  on  the 
positive  ^  axis  in  the  airfoil  plane. 

Hence  the  drag 

D  =  ?(OApdx=  -pu^  V(Ou(0  d§  . 

0  0 


[14] 
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As  shown  by  Tulin  and  Burkart  (3),  v(4)  can  be  expressed  In  terms 
of  u(4)  In  the  following  manner;  The  downwash  at  s  point  ^  on 
the  hydrofoil.  Induced  by  a  vortex  of  strength  y,  located  at  4'  is 

27r'( )  ’  Since  y  =  -2u(4'  ),  (negative  because  u(4  )  is  the 

velocity  perturbation  on  the  bottom  of  the  hydrofoil),  the  down- 
wa^h.  Integrated  over  the  ohord,  becomes 


=  v(4)  =  7  where  0  <  4  <  1 

Combining  equations  [l41  and  tl5]  and  normalizing  we  get 


[15] 


d4d4 


=  2  /^  /"  -  |(i) 
0  0  ^ 


ii+^ni-V)  * 


u(4)u(4' )d4d4 


Aair  .1 
o  J  J 


till 


4+a 


r  0  0 

and  finally,  using  equation  [3a], 

T  2 


4'+a 


d4d4' 


(a;^0) 


2 

TT 


becomes 


(a/0)  [16] 

In  the  limiting  case  of  Infinite  depth  the  drag  coefficient 

(h=oo)  [l6a] 


C  =  - 
D  TT 


u(4)d4 


L  0 


At  zero  depth  the  hydrofoil  becomes  a  planing  surface  and 
Its  cavity  drag  Is,  In  all  respects,  synonymous  with  the  spray  drag 
experienced  by  such  surfaces.  Wagner  (8)  long  ago  showed  that, 
except  In  the  immediate  neighborhood  of  the  leading  edge,  the  high 
speed  flow  past  a  planing  surface  was  Identical  to  the  flow  past 
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the  lower  surface  of  an  Identically  formed  fully  wetted  airfoil. 
The  difference  between  the  two  flows  (in  the  Immediate  neighbor¬ 
hood  of  the  leading  edge)  consists  In  the  fact  that  the  leading 
edge  suction  force  experienced  by  a  fully  wetted  airfoil  does  not 
exist  on  a  planing  surface  (or,  for  that  matter,  on  any  super- 
cavltatlng  hydrofoil)  and  Is  replaced  —  on  the  planing  surface 
only  —  by  spray.  Wagner  showed  that  the  total  drag  on  a  planing 
surface  was  equal  to  half  the  drag  on  the  Identical  fully  wetted 
airfoil  less  half  the  leading  edge  suction  force  S  acting  on  the 
latter.  Since  the  total  drag  on  a  fully  wetted  airfoil  in  an  In- 
vlscld  fluid  Is  zero,  the  spray  drag  can  be  written  as 
D  =  -  I  S  (h  =  0) 

The  magnitude  of  the  leading  edge  suction  force  S  Is  well  known 
in  airfoil  theory  and  Is  given  by  Jones  (9)  as 


S 


1  Im 

x-^  0 


-  T  p  U 

r  a 


Combining  the  last  two  equations,  normalizing  with  respect  to 
2 

and  noting  that  at  h-0,  x=^  we  get,  for  the  spray  drag  co 
efficient 


C 


D 


(h=0) 


[I6b] 


Equation  [l6b]  shows  that  If  the  pressxire  (proportional  to 
u)  on  the  lower  surface,  at  the  leading  edge  of  the  foil  goes  to 
Infinity  as  VT-  the  flat  plate  Is  a  case  In  point  --  the  spray 
drag  is  indeed  a  finite  quantity.  It  also  shows,  as  Wagner  (8) 
anticipated,  that  by  proper  use  of  camber  It  Is  possible  to  produce 
lift  without  spray  or  spray  drag.  It  will  be  seen  below  that  the 
two-term,  three-term,  five-term  and  constant  pressure  camber  con¬ 
figurations  fulfill  this  ambition. 

The  ordinates  of  the  lower  and  upper  cavity  streamlines 
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(y  and  y  respectively)  are  obtained  by  Integrating  their  slopes, 
o  c 

The  slope  of  the  lower  cavity  boundary  (l.e.,  the  bottom  surface 
of  the  foil)  has  already  been  given  In  equation  [15]*  Employing 
the  same  reasoning  but  keeping  In  mind  that  the  upper  part  of  the 
hydrofoil  Is  mapped  Into  the  negative  ^  axis,  the  slope  of  the  upper 
cavity  streamline  can  be  shown  to  be* 

^  =  V(-0  =  i  dV  where  -»  <  4  <  0  [17l 

and  the  ordinates  of  the  upper  cavity  boundary  of  the  hydrofoil  and 
lower  surface  are 

7^(0  =  V  (-!')  where  -a  <  4  <  0  [l8] 

and 

yji)  =  V  (4'  )  -■■■Al.L  av  where  0<4<1  [19] 

Rote  that  [l8]  and  [19]  give  the  cavity  ordinates  In  the  hydro¬ 
foil  plane  In  terms  of  abclssas  In  the  airfoil  plane.  This  is  only 
due  to  the  fact  that  the  transformation  (equation  [2])  cannot  be 
Inverted  to  give  an  explicit  expression  for  4  terms  of  x,  and 
can  be  remedied  In  the  process  of  plotting  the  actual  cavity  shapes. 

In  the  limiting  cases  of  zero  and  Infinite  depth  the  upper 
and  lower  cavity  ordinates  become 


* 

It  Is  pointed  out  that  in  spite  of  the  similarity  of  the  Integrands 
In  equations  [15]  and  [I7]  the  expressions  (after  Integrations)  for 
v(-4)  and  v(4)  will  be  entirely  dissimilar  due  to  the  fact  that 
equation  [15]  contains  a  singular  Integral. 
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(«')  d?' 

0 

yc(^)  =  I 

2  v(e  )  vdi 

0 

0 

2  v(?')  Vdi 

(h  =  0) 


(h  =  oo) 


[20] 


In  equations  [11],  [12],  [l6],  [l8],  and  [19]  we  have  the 
means  to  calculate  the  lift,  moment,  drag  and  upper  and  lower  cav¬ 
ity  ordinates  of  a  supercavitating  hydrofoil,  designed  for  any  spe¬ 
cified  depth,  in  terms  of  the  u  perturbation  velocity  (and  hence, 
the  pressure)  on  the  bottom  of  Its  equivalent  airfoil.  The  condi¬ 
tions  relevant  to  the  validity  of  these  equations  are  that  (a)  a 
cavity  does  Indeed  exist  eversrwhere  on  the  upper  hydrofoil  surface, 
(b  )  there  Is  no  cavitation  anywhere  on  the  lower  surface,  (c)  the 
linear  theory  holds,  (d)  the  cavitation  nxunber  is  zero,  and  (e)  the 
Proude  Nxiraber  is  very  large. 

It  is  often  convenient  to  express  the  pressures  on  a  given 
hydrofoil  as  the  sum  of  a  number  of  different  types  of  pressure 
distributions.  Conversely,  analysis  of  the  characteristics  of 
various  types  of  press\ire  distribution  often  leads  to  the  conclu¬ 
sion  that  it  might  be  advisable  to  combine  some  of  them  into  one 
foil.  Such  a  foil  is  referred  to  as  a  composite  foil  whose  various 
components  are  the  shapes  associated  with  each  of  the  specific 
types  of  pressiore  distribution.  Each  of  the  characteristics  of 
the  composite  foil  is  a  function  of  all  its  component  pressure 
distributions  in  addition  to  being  a  function  of  depth.  Denoting 
n  different  component  pressure  distributions  by  u^,  u^,  ...,  u^. 


the  characteristics  of  a  composite  foil  are  found  from  equations 
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[11],  [12],  [l6],  [l8]  and  [19]: 

^2'*  •  •  *  ~  •  •  •  ■»  ^ )  (h,  0,  ^g>  ...  0) 


+  . . .  +  Cj^(h,0,0, . .  .^u^) 


[21] 


Cj^(h,u,U2,  •  *  •-%)“C„(h,u^,0, . .  .0)+Cj^(h,0,U2,  .  .  .,0) 

■t  ...  Cjj  ( h,  0,  0,  .  . .  ) 


[22] 


C^(h,u^,u^,  . .  .  ,U^)=(  y  Cjj(h,u^,  0,  .  . . ,  0)  +yajj(Ti,0,U2,  ...,0) 

2 

+  . . .  +  -^^'c^^(it7o7^^~7TTu^^^  )  [  23  ] 

y^(h,x,u^,U2,  •  '  ••>0)+y^(h,x,0,U2,  *  ’ 

+  . . .+  y^(h,x,0,0,  ,  .  .,u^)  [24] 

and 

y^(h,x,Ui,U2, . .  .,u^)=y^(h,x,u^,0,  .  .  .,0)+y^(h,x,0,U2,  •  *  *^0) 

+  ...  +  y^(h,x,0,0,  . .  .,u  )  [25] 

^  n 

Notation  of  the  type  Cj^(h,u^,U2, .  .  .  ,u^),  ^^^roduced  here  In 

general  terms,  will  be  used  later  to  denote  the  characteristics  of 
depth-adapted  hydrofoils  which  are  combinations  of  camber  (denoted 
by  kj^),  angle  of  attack  6  and  quasi-parabolic  thickness  t.  The 

types  of  camber  configurations  that  will  be  considered  are  the 
constant  pressure  (k^),  the  two-term  (k^),  the  three-term  (k^)  and 

the  five- term  (k^).  Thus,  for  example,  y^(h,x,k2,  6,  t)  denotes 
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• 

the  bottom  surface  ordinate,  at  x,  of  a  three-term  foil  which  has 
been  designed  for  operation  at  a  depth  of  h  chords  and  whose  camber 
Index,  angle  of  attack  and  thickness  coefficient  have  the  values 
k^,  6  and  t,  respectively;  Cj^(l, 0,2.5°)  Is  the  lift  coefficient  de¬ 
veloped  by  a  flat  plate  operating  at  a  depth  of  1  chord  at  an 
angle  of  attack  of  2.5  degrees;  and  Gjj(l,0, 2,5°,  0. 01 )  Is  the  drfiig 
coefficient  of  a  quasl-parabollc  shape  designed  for  operation  at 
a  depth  of  one  chord  with  a  thickness  coefficient  of  at  an  angle 
of  attack  of  2.5  degrees. 

The  general  expressions  will  now  be  used  to  calculate  the 
lift,  drag,  moment  and  shape  of  depth  adapted  camber  configurations 
which  are  characterized  by  certain  specified  tj-pes  of  pressure 
distributions . 


THE  OPTIMUM  PRESSURE  DISTRIBUTION 

Tulin  and  Burkart  (3)  have  shown  that  for  the  class  of  all 
possible  cavltatlng  hydrofoils  the  absolute  optimum  pressure  dis¬ 
tribution  Is  that  In  which  the  pressure  Is  entirely  concentrated 
at  the  trailing  edge  of  the  foil.  They  pointed  out  that  In  spite 
of  the  fact  that  such  a  pressure  distribution  cannot  be  realized 
on  any  practical  hydrofoil  the  performance  of  this  "optimum  foil" 

Is  a  very  convenient  criterion  against  which  the  lift-drag 
efficiency  of  all  other  supercavltatlng  foils  can  be  measured.  The 
lift-drag  ratio  of  the  optimum  foil,  designed  for  infinite  depth. 

Sir 

was  shown  to  be  t;—  .  It  is  of  Interest  to  find  the  value  of  this 

optimum  ratio  for  foils  designed  for  other  depths. 

Let  us  consider  the  pressure  distribution  in  the  airfoil 
plane,  shown  In  Figure  4,  which  is  such  that  the  pressiires  are 
entirely  concentrated  at  a  point  The  distribution  of  the  H" 

velocity  on  the  bottom  s\irface  of  the  foil  can  then  be  expressed  as 
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u(e)  =  -K6(e-4o) 

where  K  is  a  constant  and  6(^-4^)  is  the  Dirac  delta  function 
which  is  such  that 


6(5-5  )  *  °  ^  ° 

°  00  for  =  0 


r  6(e-?^)d^  =  1 
.00 

and 

-.00 

Combining  equation  [26]  with  (N)  and  [l6]  we  obtain 


[27] 


Cleanly,  the  largest  possible  lift-drag  ratio  will  be  obtained 
when  4q  =  1  (l.e.,  the  pressures  are  concentrated  at  the  trailing 
edge )  so  that 


(-1  =  4h 

2_ 

4-  ^ 

opt  L 

Figure  5  Illustrates  how  the  optimum  lift-drag  ratio  for  a  given 

lift  coefficient  varies  with  depth.  As  can  be  seen,  for  a 

'^/opt 

depth  of  one  chord  is  almost  four  times  as  high  as  that  for  infinite 
depth.  Equation  [27]  also  Illustrates  that  in  the  absence  of 
strength  requirements  (which,  in  a  practical  case,  might  be  over¬ 
whelming)  the  center  of  pressure  should  be  located  as  far  aft  as 
possible . 
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THE  TWO-,  THREE-  AND  PIVE-TERM  CAMBER  DISTRIBUTIONS 

The  perturbation  velocity  In  the  airfoil  plane  on  the  two- 
terra  camber  distribution  Is  given  by  Tulin  and  Burkart  (3)  as 

u(4)  =  -  ^  k2(slne-i  sin  29)  [28] 

where  6  Is  a  dummy  variable  defined  by 
4  =  i(l-cos0) 

Johnson  (4)  showed  that  for  the  three-term  camber 

^ii)  "  ■  3^  k2(sln  9  -  sin  29  +  ^  sin  36)  [29] 

and  for  the  five-term  camber 

u(4)  =  -  k^(sln  0  -  j  sin  20  +  j  sin  30  -  j  sin  40 

+  j  sin  56)  [30] 

the  u  distributions  for  given  values  of  k^,  k_,  and  k^  are  shown 

23  5 

In  Figure  6.  The  coefficients  k^,  k_,  and  k  serve  the  purpose  of 

scaling  the  magnitude  of  the  respective  pressure  distributions  and, 

hence,  serve  also  as  scaling  factors  for  the  C^,  C„,  C_  and  shape 

L  n  D 

of  the  three  types  of  camber  configuration.  These  coefficients 

shall  be  termed  the  2-,  3-  and  5-term  camber  Index.  It  will  be 

seen  that  for  camber  configurations  designed  for  Infinite  depth 

the  numerical  value  of  the  camber  Index  Is  also  equal  to  the  value 

of  the  lift  coefficient  when  no  angle  of  attack  Is  present. 

When  the  expressions  for  u(^)  given  In  equations  [28],  [29] 

and  [30]  are  combined  with  equations  [ll],  [16],  [I8],  and  [19]  the 

following  expressions  for  the  lift,  drag,  bottom  shape  and  top 

cavity  shape  of  the  three  camber  configurations  are  obtained: 

*  ~  ■  ^  ^  “ 

The  work  Involved  Is  one  of  straightforward.  If  sometimes  laborious. 
Integrations  demanding  special  care  when  Improper  Integrals  are  In¬ 
dicated.  For  some  of  the  more  complicated  Integrations  the  tables 
of  Integrals  In  [10]  and  [11]  were  used. 
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Plrst.  let 


T  a+i  air 
In  — ^  =  •: — 


A  =  -  a  in  — ^  ~  X 

^  ah 


B  =  -  aA 

aB 


c.i?3 


D  =  i  C  -  aC 
E  =  -  aD 
F  =  I  -  al 


R 


=  i  ror  U  <  0) 


G  =  1  (4  -  i)  R  -  In  (1-2I-2R) 


H  =  ~  R^  +  I  G 


or  =  t  I  r3  +  I  H 


K  =  5  +  10 

L  ■=  +  i  it 


M  •  R  -  (a+i)ln(l-25-2R)-  "^a^+a  In  ° - ''/a^+a _R 
where  ^  Is  defined  by  equation  [2], 


let  Q  = 


g'  =  -  ^  (V^  +  i)  Q  -  I  In  (1  +  2YF  -  2Q) 


H'  =  j  +  i  g' 
j'  =  -  i  Vx  f  h' 
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M=  yxQ  +YI^L 
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and  let  b  =  2a  +  2 


+  1 


Then,  for  the  two-term  camber i 

5  2 


(h^kg^O)  =  { 


C^(h,k^,0,0)  = 


4hk, 

_ £ 

5a7r 


k. 


0 


1  +  8a^  (-yfr  -  -\fl+a  -  2a 


2 

5ir) 


25r  2 


l-4a  ■  Vi+®  )* 


(h  =  0) 

(0  <  h  <  ») 

(h  =  oo) 

(h  =  0) 

(0  <  h  <  ») 

(h  =  •) 


^^^2  ,  ,  _  2  8  3x 
5;-  (X  +  2X  -  -  X  ) 


4hk^  r 


y^(h,  X,  \z^,0,0)=( 


3t 


(l-4a-8a^)(|-  -  In 

+  2a(|-)2(l+2a  -  ^  i) 


4k 


2  ,  ^  8  3/2  ,,  2x 

(x  +  j  X  /  -  4  X  ) 


(h  =  0) 


(0  <  h  <  ») 


(h  =  ») 


where  ^  is  defined  by  equation  [2] 
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/ 


8hk 


^  (i  A  H-  2B  -  4C  -  4H  -I-  4aa  -  4a^ 


yc(h,x,k  ,0,0)  -/ 


5a  ir' 
8k 


5/  (ix-ix3/2-2x=-8  J’) 


For  the  three -term  camberi 


^  ^k 
3  ^3 


C^(h,k2,0) 


2hk_ 

L  — 2. 

'  3ajr 

k^ 


.  Is.  (1  .  i  +  -i_) 

h  ^  0/ 


2b  ‘ 


Cjj  (h,k2,0,0)  =/ 


0 

4k 

h  (  — 2- 

'  3bT 


^  97r  ^3 


f  (1 


■h  ^  2  ' 


2b  ‘ 


(  2k, 

37r 


i  (x  -  3X=  ..  x3  -  8  x"*) 


y^(h,x,k2,  0,0)=  { 


4hk. 


3aTr 


2k. 


f  (i  A  -  3B  +  -  i6d) 


(5x  -  20  +  80  -  64  x^/^) 


M)  (0<h<-) 

(h  =  00) 

(h  =  0) 

(0<h<eo) 
(h  =  00) 

(h  =  0) 

(0<h<oo) 
(h  =  o») 

(h  =  0) 

(0<h<oo) 

(h  =  00) 
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y^(h,x,k2,o,o)  =/ 


2hk, 


3ar 


A-6Bf32(C-D-J)  +  (32a+l6)  H 


+  (32a  +l6a+2)  (aM-O) 


^  1 1  +  2x^/^+8x^+  x^/^-32K' 

+  16j'-2H' 


V 


For  the  five-term  camber: 


^^k 

5  5 


C^(h,k^,0) 


SaTT 


1.  1  +  iL 

^  b^ 


b^  b 


(  0 


Cjj(h,k^,0,0) 


2k. 


<  h  (3-i  +  ^ 


^k? 


— + — )® 
b3  b** 


k  50ir  5 


t  4k 

^  {x 
dir  ' 


y^(h,x,k^,0,0) 


8x^  +  40x^ 


128  6  k 

—  X  ) 


hk 


{  (4A-64&f480C-l472D+2048E-1024P) 

5aT 

4k 

^  (^.  |2  xV2+60x2-  I2i  :,5/2^  ^3 

5ir  '  3  5  3 

-  ^  x7/2) 


(0<h<«»>) 

(h  =  «) 

(h  =  0) 

(0<h<») 

(h  =  oo) 

(h  =  0) 
(0<h<«) 

(h  =  «) 

(h  =  0) 

(0<h<») 


(h  =  «) 
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y^(h,x,k^,0,0)=( 


i  r  -  8b  +  60C  -  i84D  +  25^1  -  128(P  +  L) 
+  (l28aVl92a^+104a^+24a+2)(aM  -  O) 
(l28a3+192a2+l04a+24)  H 


(0<h<o.) 


-  (128a 2+1 92a+104 )j+(l28a+192 )K 

i  (1  ,  ^  I  x3/2+i5x"^  ifi  x5/2+  If  x3+  M  xV2 

-  128M«-2H'+24J'-104k'+192L')  (h.  = 


THE  CONSTANT- PRESSURE  CAMBER  DISTRIBUTION 
Since  the  perturbation  velocity  u  is  proportional  to  the  pres¬ 
sure  (equation  [9]),  the  u(^)  associated  with  pressures  that  are 
constant  over  the  chord  01*  the  hydro^foll  (and,  herree,  the  alrroll ) 
may  be  written 

u(^)  =  constant  ~  ~  "2 

and  Is  shovm  in  Figure  6.  When  this  u(^)  Is  combined  with  equations 
[11],  [16],  [l8]j  and  [19]  the  following  expressions  for  lift,  drag, 
bottom  shape  and  top  cavity  shape  are  obtained: 

Cj^(h,k^,0)  =  k^  (0  <  h  <  «) 


C^(h,k,,0,0)  =<h  (  ^  In  ^ 


a 


(  h  =  0) 

(0  <  h  <  ») 


(h  =  00) 


Due  to  the  standard  (NACA)  practice  of  referring  to  a  constant- 
pressure  mean  line  as  an  "a  =  1  mean  line"  k^  has  been  adopted  to 

denote  the  camber  Index  associated  with  this  type  of  pressxire  dis¬ 
tribution.  The  numerical  value  of  k^  is  the  same  as  the  value  of 

the  lift  coefficient  developed  by  this  type  of  camber  (without  angle 
of  attack)  at  any  depth. 
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[x  In  X  +  (l-x)  In  (l-x)  ] 


(h  =  0) 


y^(h^x,  k^,0,0)= 


2Tr' 


(l+a)2 


_1 

2Tr 


(x-l)  In  (l--^Gc)  -  X  In  ~ 


(0<h<eo) 
(h  =  oo) 


y^(h,x,k^,0,0)= 


hk,  r 


— I  ll  in  ^  i  In  (1-0  -  In  ^  In  (l+a) 


27r 


2  La  ^ 


1 

P-TT 


(x-l)  In  (l+Y^)  -  X  In  "(/x  +  X 


(0<h<») 
(h  =  ») 


»  n 


w 

where  (w)  =  dllogarlthm  of  w  =  ^  for 


ri=l  n 


w 


<  1 


F" 


and  ^  is  defined  by  equation  [2]. 

Due  to  the  inordinate  difficulties  Involved  in  integrating 
equation  [12]  analytically,  the  pitching  moments  were  calculated 
by  performing  a  separate  numerical  integration  for  each  depth  and 
each  type  of  pressure  distribution.  This  was  accomplished  through 
the  use  of  the  IBM  1620  digital  computer.  Installed  at  HYDRONAUTICS, 
Incorporated. 

The  variation  of  lift  coefficient,  drag  coefficient  and  the 
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location  of  the  center  of  pressure  (  77-  )  with  depth,  for  the  two-, 

three-,  and  five-term  and  the  constant  pressure  types  of  camber  are 
shown  in  Figures  7(a),  7(b),  and  7(c),  respectively.  As  a  typical 
example,  the  variation  of  the  lift-drag  ratio  of  the  twxi- term  cam¬ 
ber  with  lift  coefficient,  is  shown  in  Figure  8.  The  bottom  and 
top  cavity  shapes  of  the  four  types  of  camber  distribution,  for 
depth-chord  ratios  of  1  and  infinity  are  shown  in  Figui’es  9(a)  and 
9(b),  respectively.  Except  in  the  case  of  the  constant  pressure 

type  of  camber  distribution,  the  upper  cavity  boundary  y  is  seen 

c 

to  lie  below  y^,  the  lower  -  or  pressure  -  surface  of  the  foil. 

This  result,  which  implies  a  crossing  of  the  streamlines,  is  not 
attributed  to  inadequacies  of  the  linear  theory.  The  same  pheno¬ 
menon  occurs  in  many  cases  in  which  body  shapes  are  derived,  as 
is  done  here,  from  prescribed  pressure  distributions,  regardless 
of  whether  the  problem  is  treated  with  linear  or  exact  theory  (see 
(12)  and  (13)  for  discussions  of  this  subject  in  the  literature). 
The  prescribed  pressure  distribution,  which  is  the  starting  point 
of  the  foil  design,  immediately  implies  a  theoretical  foil  shape 
but  the  theory,  be  it  linear  or  exact,  which  is  used  to  calculate 
that  shape,  does  not  necessarily  guarantee  it  to  be  physically 
feasible.  The  theory  does,  in  fact,  penalize  us  for  using  overly 
efficient  camber  distributions  (see  Figure  8  )  by  compelling  us  to 
accept  the  reductions  in  lift-drag  ratio  that  occur'  when  angle  of 
attack  and/or  quasi-parabolic  thickness  are  combined  with  camber 
in  order  to  make  the  foil  a  practical  possibility  and  provide  it 
with  sufficient  strength. 


4 
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QUASI-PARABOLIC  THICKNESS 

A  vortex  of  strength  y  =  ir  T  located  at  the  leading  edge 

of  the  equivalent  airfoil  will  generate  a  pressure  concentrated  en¬ 
tirely  at  the  leading  edge  of  the  hydrofoil  so  that  the  associated 
u  -velocity  distribution  may  be  written  In  terms  of  the  Dlrac-delta 
function  (equation  [26])  as 


U(0  -  -  I  T  6  {«) 

Conrblnlng  u(?)  with  equations  [11],  [l6],  [l8],  and  [19] 
yields  the  following  expressions  for  the  lift,  drag,  moment  and  shape 
associated  with  the  leading-edge-concentrated  pressure; 

(0  <  h  ^  «) 

(0  1  h  1  ») 


C  due  to  T  =0 
L> 


C„  due  to  T  =0 
M 


Cj^(h,0,  0,t) 


TT  2 
2 


(0  <  h  <  oo) 
(h  =  oo) 


^  (0,0,0,t)  =  0 


(h  =  0) 


y^(h,x,0,0,T)  = 


y^(h,x,  0,0,t) 


h 

2Tra 

In 

j|j± 

a 

(0 

< 

h  < 

QO) 

T  Yx 

(h 

= 

oo) 

h 

^  27ra 

In 

l+a 

a 

(0 

< 

h  < 

oo) 

(h 

= 

oo) 

[2]. 

The  effect 

of 

depth 

on 
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Is  shown  In  Figure  7(b).  The  upper  and  lower  cavity  shapes  due 
to  T,  for  a  number  of  depth-chord  ratios,  are  shown  In  Figure  10. 

THE  INVERSE  PROBLEM 

To  find  the  characteristics  of  a  supercavltatlng  hydrofoil 
of  given  shape,  using  the  methnds  outlined  herein,  it  Is  necessary 
to  express  the  pressure  distribution  -  or  u(4)  -  In  terms  of  the 
shape  before  the  expressions  (equations  [ll],  [12],  and  [16])  for 
the  lift,  moment  and  drag  coefficients  can  be  used.  The  Inversion 
of  the  integral  equation  In  [15]  was  shown  by  Munk  (7)  to  be 

u(e)  =  -  —  dc',  0  ^  e  <  1  [311 

ttVT  0  Vl-^'  (4  -  ■ 

In  which  v(4')  can  be  found  from  the  shape  by  using  equation  [13]. 

Unfortunately  the  Integrations  Involved  In  these  operations 
when  the  hydrofoil  Is  operating  at  a  finite  depth  are,  for  all  but 
the  very  simplest  shapes,  too  complicated  to  be  done  analytically 
and  would  have  to  be  done  by  numerical  methods.  At  present,  only 
the  flat  plate  problem  will  be  solved  for  finite  depths.  The  prob¬ 
lem  of  a  flat  plate  with  flap  near  a  free  surface  has  been  dealt 
with  elsewhere  (l4)  and  for  very  close  approximations  to  the  per¬ 
formance  characteristics  of  two-,  three-,  and  five-term  and  circular 
arc  foils  which  are  designed  for  infinite  depth  but  operating  at 
some  other  depth,  the  reader  Is  referred  to  Johnson's  work  (6). 

The  Flat  Plate  at  Finite  Depths 

flat  plate  at  an  angle  of  attack  6  to  the 
Is 


The  slope  of  a 
direction  of  the  flow 


I 
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We  use  equations  [13]  and  [31]  to  obtain 


.vr  -b  ^  (e- 


eo  VT 


[32] 


and  combine  this  result  with  equations  [11],  [I6]  and  [19]  to  get, 
for  a  flat  plate, 

/ 

6t  (h  =  0 

Vi  2 

(h,0,6)  =  /  6  -  (1/  1+a  -  ^JT  )  (0  <  h  <  «) 

I  \  a 

(h  =  00) 


'I 


'^7r6^ 


Cjj(h,  0,6,0)  =  / 


8^  §  (VT+a  -  Ya”  )“  (0  <  h  <  00) 


E 

2  ® 


(h  «  0) 


(h  =  co) 


^6^ 
TT 


y^(h,x,0,6,0)={ 


-  2x 


(0<h<<«>) 


where 


{l+2^fx  )-y^x+Y3r 

^  In  (1+  2Yx  -  2-^yi+^fT  ) 

P  =  In  (1  -  2e  -  2-yPT ) 


(h  =  «) 


i+& 


Q  =  in  (1  -  ^  §  +  2y^  )  -  In  g 

and  I  is  defined  by  equation  [2]. 
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The  variation  of  Cj^(h,  0, 6)  and  Cp(h,  0, 6>0)  and  Cj|j(h,  0, 6) 

with  depth  Is  shown  In  Figures  7(s)  7(t).  The  moment  coef¬ 

ficients  were  calculated  by  combining  equations  [32]  and  [11]  and 
performing  a  numerical  Integration.  The  result  Is  shovm  In  Fig¬ 
ure  7(c)  as  a  plot  of  location  of  the  center  of  pressure 

C|y. 

(■~)  on  the  hydrofoil  vs.  depth-chord  ratio.  The  top  cavity  shape 

due  to  angle  of  attack  for  a  number  of  depth-chord  ratios  Is  shown 
In  Figure  11. 

OPTIMUM  FOILS  AND  FOIL  STRENGTH 

In  accordance  with  equations  [21],  [23],  [24],  and  [25]  the 
lift,  drag  and  upper  and  lower  cavity  ordinates  of  a  foil  In  which 
camber,  angle  of  attack  and  quasi-parobollc  thickness  are  combined 
are  found  as  follows: 

Cj^(h,k^,6)  =  Cj^(h,k^,0,  )  +  C^(h,0,6) 

CD(h,ki,6,T)  =  |yCjj(h,k^,0,0)  +  ^  C^{h,0,b,0)  + 

Cj^(h,k^,6)  =  Cjj(h,k^,0)  +  Cjj(h,0,6) 

y^(h,x,k^,6,T)  =  y^(h,x,k^,0,  0)+y^(h,x,  0,6,0 )+y^(h,x,  0,0, t) 

y^(h,x,k^,6,T)  =  y^(h,x,k^,0,0)+y^{h,x,0,6,0)+y^(h,x,0,0,T) 

Of  the  four  camber  configurations  considered  here  It  is 
clear  that  for  given  values  of  6  and  x,  foils  Incorporating  the 
five-term  camber  produce  the  highest  lift-drag  ratios.  This  con¬ 
clusion  is  no  longer  significant,  however,  when  the  strength  of 
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these  supercavltatlng  hydrofoils  Is  taken  Into  consideration.  For 
purposes  of  strength  calculations  it  is  assumed  that  the  foil  sec¬ 
tion  is  constructed  so  as  to  fill  the  entire  thickness  of  the 
cavity*,  thus  setting  the  foil  thickness  at  any  point  along  the 
chord  equal  to  y^(h,x,  k^,  6,t)  -  y^(h,x, k^,  6,t ) .  The  foil  strength 

is  evaluated  by  considering  the  non-dimensional  section  modulus  Z 
of  a  beam  whose  cross  section  is  the  hydrofoil  section  [z  =  Z/(c)^ 
where  Z  =  (area  moment  of  inertia) ^(the  distance  between  the  neutral 
axis  and  the  furthest  flber)j.  In  these  calculations  it  is  assumed 
that  the  neutral  axis  of  the  beam  is  parallel  to  the  reference 
line  (y  =  O)  of  the  foil.  Since  the  section  modulus  is  calculated 
mainly  for  comparison  purposes  it  is  of  little  consequence  that 
the  true  neutral  axis  may  be  oriented  somewhat  differently  and  that, 
in  addition,  it  may  not  always  be  practical  to  make  the  foil  quite 
as  thick  as  the  cavity. 

The  integrations  Involved  in  calculating  the  section  modulus 
must,  of  necessity,  be  carried  out  numerically  for  each  specific 
foil  (l.e.,  a  foil  composed  of  a  given  set  of  the  parameters  h,k^, 

6  and  t).  These  calculations  were  first  carried  out  by  hand  for 
the  two-term  foil  at  a  depth  of  one  chord,  using  Simpson's  Rule 
in  the  calculation  of  the  section  moduli.  Subsequently,  a  high 
speed  digital  computer  at  the  David  Taylor  Model  Basin  was  used  to 
calculate  the  operating  characteristics  and  section  moduli  of  two- 
term,  three-term,  five-term,  and  constant  pressure  camber  hydrofoils, 
designed  for  depth-chord  ratios  of  1/4,  1/2,  1,  2,  5  and  ».  in 
these  calculations  the  computer  was  programmed  to  consider  for  each 
depth  (a)  all  possible  combinations  of  k^^,  6,  and  t,  with  k^^^  taking 


*The  upper  surface  of  the  foil  is,  however,  still  considered  to  be 
unwetted. 
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on  10  values  between  0  and  0.6,  6  taking  on  13  values  between  0  and 
3  degrees,  and  t  taking  on  9  values  between  0  and  1^;  and  (b)  all 
possible  combinations  of  the  constant  pressure  camber  configurations 
with  each  of  the  two-,  three-,  and  five-term  camber  configurations, 
with  k^  taking  on  the  same  values  as  In  (a).  A  sample  of  the  re¬ 
sults  of  these  calculations  Is  shown  In  Figure  12,  In  which  con¬ 
tours  of  constant  6,  k.^  and  Z  for  a  single  value  of  t  are  plotted 
on  a  grid  of  lift-drag  ratio  vs.  lift  coefficient,  for  a  two-term 
hydrofoil  designed  for  a  depth  of  one  chord.  A  subsequent  report 
(No.  001-7)  containing  all  the  significant  results  of  this  computer 
study  will  be  published  In  the  near  future.  In  addition,  the  com¬ 
puter  program  contained  an  optimization  procedure  In  which  the  pre¬ 
viously  stored  list  of  all  the  foils  that  had  been  computed  for  a 
given  depth  was  searched  for  the  three  best  foils  (l.e,,  those 
with  highest  L/D)  whose  section  modulus  and  lift  coefficient  were 
within  certain  prescribed  ranges. 

The  optimization  procedure  showed  that  It  Is  not  possible 
to  single  out  one  type  of  camber  distribution  as  being  the  absolute 
optimum  for  all  depths,  lift  coefficients  and  section  moduli.  It 
Is  possible,  however,  to  make  the  following  generalizations: 

(a)  The  combinations  of  the  constant  pressure  camber  with 
two-,  three-,  and  five-term  camber  configurations  are  to  be  recom¬ 
mended  only  at  very  low  depth-chord  ratios  (h  <  0.5).  At  depth- 

chord  ratios  of  0.5  or  higher  this  type  of  foil  Is  efficient  only 

—  -4 

for  unrealistically  high  section  moduli  (Z  >  12  x  10  )  and  un¬ 

realistically  high  lift  coefficients  (C  >  .60). 

ij 

(b)  Turning  now  to  foils  In  which  one  of  the  four  camber 
types  Is  combined  with  angle  of  attack  and  quasl-parabollc  thick¬ 
ness,  It  appears  that  except  at  unrealistically  low  section  moduli 
(Z  <  2  X  10  ^)  the  constant  pressure  and  two- term  are  more  efficient 
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than  the  three-  and  five-term  camber  configurations.  This  trend 

Intensifies  as  either  the  section  modulus  oi^  the  design  depth  Is 

increased,  until,  at  infinite  depth  and  high  section  moduli 

(Z  >  6  X  10  the  constant  pressure  type  of  camber  (combined  with 

6  and  t)  dominates  the  list  of  optimum  foils. 

In  summary,  the  optimization  process  Indicates . that  for 

practical  values  of  depth-chord  ratio,  section  modulus  and  lift 

_]i  —  -4 

coefficient  (i.e. ,  h  >  0.5>  2  x  10  <  Z  <  6  x  10  and  0.05  < 

C^.  <  0.5)  the  constant  pressure  and  two-term  camber  configurations, 
each  properly  combined-  with  angle  of  attack  and  quasi-parabolic 
thickness,  deliver  higher  lift-drag  ratios  than  the  other  types  of 
camber.  The  two- term  foils  are,  usually,  the  better  of  these  two 
and  in  cases  In  which  the  constant  pressure  camber  Is  somewhat 
superior  the  difference  between  Its  3ift-drag  ratio  and  that  of 
the  two-term  foils  is  not  very  significant. 

In  Figure  13  a  typical,  two-dimensional,  two-term  hydrofoil 
is  drawn  to  scale.  This  particular  foil  is  designed  for  operation 
at  a  depth  of  one  chord  and  is  composed  of  camber  Index  k^  =  O.15, 
design  angle  of  attack  B  =  1.6°  and  quasi-parabolic  thickness 
T  -  0.004.  It  develops  a  lift  coefficient  of  0.195  and  a  lift- 

cavity  drag  ratio  of  34.5  and  Its  nominal  section  modulus  is 

—  -4 

Z  =  4,1  X  10  . 
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Figure  3  -  Transformation  Parameters 
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Figure  7-  The  influence  of  Depth  on  Lift  Coefficient,  Drag  Coefficient 
and  Location  of  the  Center  of  Pressure  due  to  Camber,  Angle 
of  Attack  and  Quasi -Parabolic  Thickness 
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Figure  10-  Bottom  Surface,  ,  and  Top  Cavity  Boundary,  ,  of  Quasi-  Parabolic 

Thickness,  at  Various  Depths  of  Submersion 
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Figure  12-  Lift-Drag  Ratio  versus  Lift  Coefficient  of  a  Two-Term  Foil  with 
Angle  of  Attack  and  Quasi-Parabolic  Thickness, ^"=0.004  •,h  =  I 
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